Bayesian Inference Of Phylogeny Research Articles

A New Sturgeon Herpesvirus from Juvenile Lake Sturgeon Acipenser fulvescens Displaying Epithelial Skin Lesions.

Herpesvirus infections of sturgeon pose a potential threat to sturgeon culture efforts worldwide. A new epitheliotropic herpesvirus named Acipenser herpesvirus 3 (AciHV-3) was detected in hatchery-reared Lake Sturgeon Acipenser fulvescens displaying skin lesions in central Canada. The growths were discovered in the fall, reached average prevalence levels of 0.2-40% and eventually regressed. No unusual mortality was observed. The cellular changes within the lesions included epithelial hyperplasia and were reminiscent of other herpesvirus infections. The virus was not evident in lesions examined by electron microscopy. Skin tissue homogenates from symptomatic sturgeon produced atypical cytopathic effects on a primary Lake Sturgeon cell line, and next-generation sequence analysis of the DNA samples revealed the presence of an alloherpesvirus. A new genotyping PCR assay targeting the major capsid protein sequence detected AciHV-3 in symptomatic Lake Sturgeon as well as other apparently healthy sturgeon species. Bayesian inference of phylogeny reconstructed with a concatenation of five alloherpesvirus core proteins revealed a new Alloherpesviridae lineage isomorphic with a new genus. The presence of AciHV-3 homologs in cell lines and sturgeon sequence datasets, low sequence divergence among these homologs and branching patterns within the genotyping phylogeny provide preliminary evidence of an endogenous virus lifestyle established in an ancestral sturgeon.

Nucleotide Substitution Model Selection Is Not Necessary for Bayesian Inference of Phylogeny With Well-Behaved Priors.

Model selection aims to choose the most adequate model for the statistical analysis at hand. The model must be complex enough to capture the complexity of the data but should be simple enough not to overfit. In phylogenetics, the most common model selection scenario concerns selecting an adequate substitution and partition model for sequence evolution to infer a phylogenetic tree. Previously, several studies showed that substitution model under-parameterization can bias phylogenetic studies. Here, we explored the impact of substitution model over-parameterization in a Bayesian statistical framework. We performed simulations under the simplest substitution model, the Jukes-Cantor model, and compare posterior estimates of phylogenetic tree topologies and tree length under the true model to the most complex model, the $\text{GTR}+\Gamma+\text{I}$ substitution model, including over-splitting the data into additional subsets (i.e., applying partitioned models). We explored 4 choices of prior distributions: the default substitution model priors of MrBayes, BEAST2, and RevBayes and a newly devised prior choice (Tame). Our results show that Bayesian inference of phylogeny is robust to substitution model over-parameterization and over-partitioning but only under our new prior settings. All 3 current default priors introduced biases for the estimated tree length. We conclude that substitution and partition model selection are superfluous steps in Bayesian phylogenetic inference pipelines if well-behaved prior distributions are applied and more effort should focus on more complex and biologically realistic substitution models.

Thecate stem medusozoans (Cnidaria) from the early Cambrian Chengjiang biota

AbstractCnidarians are phylogenetically located near the base of the ‘tree of animals’, and their early evolution had a profound impact on the rise of bilaterians. However, the early diversity and phylogeny of this ‘lowly’ metazoan clade has hitherto been enigmatic. Fortunately, cnidarian fossils from the early Cambrian could provide key insights into their evolutionary history. Here, based on a scrutiny of the purported hyolith Burithes yunnanensis Hou et al. from the early Cambrian Chengjiang biota in South China, we reveal that this species shows characters distinct from those typical of hyoliths, not least a funnel‐shaped gastrovascular system with a single opening, a whorl of tentacles surrounding the mouth, and the lack of an operculum. These characters suggest a great deviation from the original definition of the genus Burithes, and a closer affinity with cnidarians. We therefore reassign the material to a new genus: Palaeoconotuba. Bayesian inference of phylogeny based on new anatomical traits identifies a new clade, including Palaeoconotuba and Cambrorhytium, as a stem group of sessile medusozoan cnidarians that are united by the synapomorphies of developing an organic conical theca and a funnel‐like gastrovascular system. This study unveils a stem lineage of medusozoans that evolved a lifelong conical theca in the early Cambrian.

The evolution of unique cranial traits in leporid lagomorphs.

The leporid lagomorphs (rabbits and hares) are adapted to running and leaping (some more than others) and consequently have unique anatomical features that distinguish them from ochotonid lagomorphs (pikas) and from their rodent relatives. Two traits that have received some attention are fenestration of the lateral wall of the maxilla and facial tilting. These features are known to correlate with specialised locomotory form in that the faster running species will generally have fenestration that occupies the dorsal and the anteroventral surface of the maxillary corpus and a more acute facial tilt angle. Another feature is an intracranial joint that circumscribes the back of the skull, thought to facilitate skull mobility. This joint separates the anterior portion of the cranium (including the dentition, rostrum and orbit) from the posterior portion of the cranium (which encompasses the occipital and the auditory complex). Aside from the observation that the intracranial joint is absent in pikas (generalist locomotors) and appears more elaborate in genera with cursorial and saltatorial locomotory habits, the evolutionary history, biomechanical function and comparative anatomy of this feature in leporids lacks a comprehensive evaluation. The present work analysed the intracranial joint, facial tilting and lateral fenestration of the wall of the maxilla in the context of leporid evolutionary history using a Bayesian inference of phylogeny (18 genera, 23 species) and ancestral state reconstruction. These methods were used to gather information about the likelihood of the presence of these three traits in ancestral groups. Our phylogenetic analyses found it likely that the last common ancestor of living leporids had some facial tilting, but that the last common ancestor of all lagomorphs included in the dataset did not. We found that it was likely that the last common ancestor of living leporids had fenestration that occupies the dorsal, but not the anteroventral, surface of the maxillary corpus. We also found it likely that the last common ancestor of living leporids had an intracranial joint, but that the last common ancestor of all living lagomorphs did not. These findings provide a broader context to further studies of evolutionary history and will help inform the formulation and testing of functional hypotheses.

Bayesian inference of phylogeny is robust to substitution model over-parameterization

Bayesian inference of phylogeny is robust to substitution model over-parameterization

Bayesian inference of phylogeny is robust to substitution model over-parameterization

Bayesian inference of phylogeny is robust to substitution model over-parameterization

Crouching TIGER, hidden structure: Exploring the nature of linguistic data using TIGER values

AbstractIn recent years, techniques such as Bayesian inference of phylogeny have become a standard part of the quantitative linguistic toolkit. While these tools successfully model the tree-like component of a linguistic dataset, real-world datasets generally include a combination of tree-like and nontree-like signals. Alongside developing techniques for modeling nontree-like data, an important requirement for future quantitative work is to build a principled understanding of this structural complexity of linguistic datasets. Some techniques exist for exploring the general structure of a linguistic dataset, such as NeighborNets, δ scores, and Q-residuals; however, these methods are not without limitations or drawbacks. In general, the question of what kinds of historical structure a linguistic dataset can contain and how these might be detected or measured remains critically underexplored from an objective, quantitative perspective. In this article, we propose TIGER values, a metric that estimates the internal consistency of a genetic dataset, as an additional metric for assessing how tree-like a linguistic dataset is. We use TIGER values to explore simulated language data ranging from very tree-like to completely unstructured, and also use them to analyze a cognate-coded basic vocabulary dataset of Uralic languages. As a point of comparison for the TIGER values, we also explore the same data using δ scores, Q-residuals, and NeighborNets. Our results suggest that TIGER values are capable of both ranking tree-like datasets according to their degree of treelikeness, as well as distinguishing datasets with tree-like structure from datasets with a nontree-like structure. Consequently, we argue that TIGER values serve as a useful metric for measuring the historical heterogeneity of datasets. Our results also highlight the complexities in measuring treelikeness from linguistic data, and how the metrics approach this question from different perspectives.

The first third-generation HIV-1 circulating recombinant form (CRF114_0155) identified in central China.

A novel circulating recombinant form (CRF) was identified in eight HIV-1-infected patients without direct epidemiological relationships in Henan Province, Central China. Recombination analysis indicated that the genome of this novel CRF comprises five segments: three inherited from CRF01_AE cluster-4 and two from CRF55_01B. Therefore, the CRF was designated CRF114_0155. It is not only the first novel CRF identified in Henan Province but also the first third-generation CRF of HIV-1 and the first CRF descendant of CRF55_01B. Bayesian inference of phylogeny dated the most recent common ancestor of the CRF114_0155 cluster to 2010. The emergence of CRF114_0155 reflects that the genotype constitution of HIV-1 has become more complex and that stricter intervention measures should be implemented in central China.

Diversity of Trametes (Polyporales, Basidiomycota) in tropical Benin and description of new species Trametes parvispora.

Trametes is a globally distributed genus of white-rot polypores and well sampled in temperate and boreal areas. However, the diversity, taxonomy, and phylogenetic positions of Trametes spp. are poorly known in tropical Africa. This study aims at documenting the diversity of Trametes species in Benin (tropical Africa) and their phylogenetic positions with a focus on the T. elegans species complex. Therefore, we collected specimens of Trametes from different forest types across Benin. To infer phylogenetic relationships between Trametes species, we investigated sequences of five gene regions and added available sequences from GenBank. Using Maximum likelihood and Bayesian phylogeny inference methods, we found eight supported species clades. For the T. elegans species complex, we re-establish the name Trametes palisotii for species previously known as T. elegans in tropical Africa. Furthermore, we propose Trametes parvispora as a species new to science and provide the description of this species. Our molecular phylogeny of Trametes with a focus on tropical Benin contributes to taxonomic clarity of an important wood-decay fungal genus, which is the basis for biodiversity assessments of Trametes in the tropics.

Simultaneous Bayesian inference of phylogeny and molecular coevolution

Patterns of molecular coevolution can reveal structural and functional constraints within or among organic molecules. These patterns are better understood when considering the underlying evolutionary process, which enables us to disentangle the signal of the dependent evolution of sites (coevolution) from the effects of shared ancestry of genes. Conversely, disregarding the dependent evolution of sites when studying the history of genes negatively impacts the accuracy of the inferred phylogenetic trees. Although molecular coevolution and phylogenetic history are interdependent, analyses of the two processes are conducted separately, a choice dictated by computational convenience, but at the expense of accuracy. We present a Bayesian method and associated software to infer how many and which sites of an alignment evolve according to an independent or a pairwise dependent evolutionary process, and to simultaneously estimate the phylogenetic relationships among sequences. We validate our method on synthetic datasets and challenge our predictions of coevolution on the 16S rRNA molecule by comparing them with its known molecular structure. Finally, we assess the accuracy of phylogenetic trees inferred under the assumption of independence among sites using synthetic datasets, the 16S rRNA molecule and 10 additional alignments of protein-coding genes of eukaryotes. Our results demonstrate that inferring phylogenetic trees while accounting for dependent site evolution significantly impacts the estimates of the phylogeny and the evolutionary process.

Molecular systematics of sturgeon nucleocytoplasmic large DNA viruses

Namao virus (NV) is a sturgeon nucleocytoplasmic large DNA virus (sNCLDV) that can cause a lethal disease of the integumentary system in lake sturgeon Acipenser fulvescens. As a group, the sNCLDV have not been assigned to any currently recognized taxonomic family of viruses. In this study, a data set of NV DNA sequences was generated and assembled as two non-overlapping contigs of 306,448 bp and then used to conduct a comprehensive systematics analysis using Bayesian inference of phylogeny for NV, other sNCLDV and representative members of six families of the NCLDV superfamily. The phylogeny of NV was reconstructed using protein homologues encoded by nine nucleocytoplasmic virus orthologous genes (NCVOGs): NCVOG0022 – mcp, NCVOG0038 – DNA polymerase B elongation subunit, NCVOG0076 – VV A18-type helicase, NCVOG0249 – VV A32-type ATPase, NCVOG0262 – AL2 VLTF3-like transcription factor, NCVOG0271 – RNA polymerase II subunit II, NCVOG0274 – RNA polymerase II subunit I, NCVOG0276 – ribonucleotide reductase small subunit and NCVOG1117 – mRNA capping enzyme. The accuracy of our phylogenetic method was evaluated using a combination of Bayesian statistical analysis and congruence analysis. Stable tree topologies were obtained with data sets differing in target molecule(s), sequence length and taxa. Congruent topologies were obtained in phylogenies constructed using individual protein data sets. The major capsid protein phylogeny inferred that ten representative sNCLDV form a monophyletic group comprised of four lineages within a polyphyletic Mimi-Phycodnaviridae group of taxa. Overall, the analyses revealed that Namao virus is a member of the Mimiviridae family with strong and consistent support for a clade containing NV and CroV as sister taxa.

Posterior Summarization in Bayesian Phylogenetics Using Tracer 1.7

Bayesian inference of phylogeny using Markov chain Monte Carlo (MCMC) plays a central role in understanding evolutionary history from molecular sequence data. Visualizing and analyzing the MCMC-generated samples from the posterior distribution is a key step in any non-trivial Bayesian inference. We present the software package Tracer (version 1.7) for visualizing and analyzing the MCMC trace files generated through Bayesian phylogenetic inference. Tracer provides kernel density estimation, multivariate visualization, demographic trajectory reconstruction, conditional posterior distribution summary, and more. Tracer is open-source and available at http://beast.community/tracer.

VMCMC: a graphical and statistical analysis tool for Markov chain Monte Carlo traces

BackgroundMCMC-based methods are important for Bayesian inference of phylogeny and related parameters. Although being computationally expensive, MCMC yields estimates of posterior distributions that are useful for estimating parameter values and are easy to use in subsequent analysis. There are, however, sometimes practical difficulties with MCMC, relating to convergence assessment and determining burn-in, especially in large-scale analyses. Currently, multiple software are required to perform, e.g., convergence, mixing and interactive exploration of both continuous and tree parameters.ResultsWe have written a software called VMCMC to simplify post-processing of MCMC traces with, for example, automatic burn-in estimation. VMCMC can also be used both as a GUI-based application, supporting interactive exploration, and as a command-line tool suitable for automated pipelines.ConclusionsVMCMC is a free software available under the New BSD License. Executable jar files, tutorial manual and source code can be downloaded from https://bitbucket.org/rhali/visualmcmc/.

Identification and Characterization of Rhizoctonia Species Associated with Soybean Seedling Disease.

In an effort to identify the Rhizoctonia spp. associated with seedling diseases of soybean, Rhizoctonia isolates were recovered from soybean seedlings with damping off and root and hypocotyl rot symptoms from Arkansas, Illinois, Kansas, Michigan, Minnesota, and the Canadian province of Ontario between 2012 and 2014. Based on cultural morphology, polymerase chain reaction restriction fragment length polymorphism, and phylogenetic analysis of the internal transcribed spacer (ITS) region of the ribosomal RNA genes, 80 isolates were confirmed to be Rhizoctonia solani, 24 were binucleate Rhizoctonia spp., and 10 were R. zeae. Of the 80 R. solani isolates, one belonged to anastomosis group (AG) 2-1, 52 belonged to AG-2-2IIIB, five belonged to AG-3 PT, three belonged to AG-4 HGI, two belonged to AG-4 HGIII, nine belonged to AG-7, and eight belonged to AG-11. Bayesian inference of phylogeny using the ITS region revealed two clades of R. solani AG-7 that possibly correspond to different AG-7 subgroups. Phylogenetic analysis also provided evidence for genetic relatedness between certain binucleate Rhizoctonia and some R. solani isolates. On 'Williams 82' soybean, isolates of AG-2-2IIIB were the most aggressive, followed by isolates of AG-7, AG-4, and AG-11. On 'Jubilee', a sweet corn cultivar, AG-2-2IIIB and AG-4 isolates caused significant stunting and root damage, whereas the damage caused by the AG-11 isolates was mostly restricted to the mesocotyl. Isolates of R. zeae and the binucleate Rhizoctonia spp. were not pathogenic on soybean or corn. Our results indicate that soybean and corn are hosts to the predominant and aggressive AG of R. solani, implying that rotation between these two crops may not be an effective management practice.

Morphological comparison and DNA barcoding of four closely related species in the genera Misgurnus and Paramisgurnus (Cypriniformes: Cobitidae)

To evaluate the feasibility of morphological and genetic identification of the closely related species in the genera Misgurnus and Paramisgurnus, the morphological characters of four species in these genera and DNA barcoding of five loaches (P. dabryanus, M. anguillicaudatus, M. bipartitus, M. mohoity, and Barbatula toni) were investigated. Twelve morphological characters were measured in 542 individuals to perform the comparative analysis. Among these characters, only the caudal peduncle length (LCP) revealed significant difference (P < 0.05) among these four species. The clustering based on morphological characters formed two clusters (P. dabryanus and M. anguillicaudatus; M. bipartitus and M. mohoity). A total of 186 COI fragments for the five loaches investigated were sequenced and analyzed. The results showed that interspecific K2P distance was much higher than intraspecific distance within the five species. Bayesian inference of phylogeny showed that individuals of these species were divided into five specific clades. Meanwhile, the COI fragments exhibited 22 character attributes for the differentiation of the five loach species based on character-based method. Our results suggested that DNA barcoding based on COI can be used as an efficient identifier of these five loach species; the combination of distance-based method, Bayesian inference and character-based approach provides higher resolution of identification at species level.

Clonal Evolution of Enterocytozoon bieneusi Populations in Swine and Genetic Differentiation in Subpopulations between Isolates from Swine and Humans.

Enterocytozoon bieneusi is a widespread parasite with high genetic diversity among hosts. Its natural reservoir remains elusive and data on population structure are available only in isolates from primates. Here we describe a population genetic study of 101 E. bieneusi isolates from pigs using sequence analysis of the ribosomal internal transcribed spacer (ITS) and four mini- and microsatellite markers. The presence of strong linkage disequilibrium (LD) and limited genetic recombination indicated a clonal structure for the population. Bayesian inference of phylogeny, structural analysis, and principal coordinates analysis separated the overall population into three subpopulations (SP3 to SP5) with genetic segregation of the isolates at some geographic level. Comparative analysis showed the differentiation of SP3 to SP5 from the two known E. bieneusi subpopulations (SP1 and SP2) from primates. The placement of a human E. bieneusi isolate in pig subpopulation SP4 supported the zoonotic potential of some E. bieneusi isolates. Network analysis showed directed evolution of SP5 to SP3/SP4 and SP1 to SP2. The high LD and low number of inferred recombination events are consistent with the possibility of host adaptation in SP2, SP3, and SP4. In contrast, the reduced LD and high genetic diversity in SP1 and SP5 might be results of broad host range and adaptation to new host environment. The data provide evidence of the potential occurrence of host adaptation in some of E. bieneusi isolates that belong to the zoonotic ITS Group 1.

RevBayes: Bayesian Phylogenetic Inference Using Graphical Models and an Interactive Model-Specification Language.

Programs for Bayesian inference of phylogeny currently implement a unique and fixed suite of models. Consequently, users of these software packages are simultaneously forced to use a number of programs for a given study, while also lacking the freedom to explore models that have not been implemented by the developers of those programs. We developed a new open-source software package, RevBayes, to address these problems. RevBayes is entirely based on probabilistic graphical models, a powerful generic framework for specifying and analyzing statistical models. Phylogenetic-graphical models can be specified interactively in RevBayes, piece by piece, using a new succinct and intuitive language called Rev. Rev is similar to the R language and the BUGS model-specification language, and should be easy to learn for most users. The strength of RevBayes is the simplicity with which one can design, specify, and implement new and complex models. Fortunately, this tremendous flexibility does not come at the cost of slower computation; as we demonstrate, RevBayes outperforms competing software for several standard analyses. Compared with other programs, RevBayes has fewer black-box elements. Users need to explicitly specify each part of the model and analysis. Although this explicitness may initially be unfamiliar, we are convinced that this transparency will improve understanding of phylogenetic models in our field. Moreover, it will motivate the search for improvements to existing methods by brazenly exposing the model choices that we make to critical scrutiny. RevBayes is freely available at http://www.RevBayes.com. [Bayesian inference; Graphical models; MCMC; statistical phylogenetics.]

Clinicopathologic, Immunohistochemical, and Ultrastructural Findings of a Fatal Case of Middle East Respiratory Syndrome Coronavirus Infection in the United Arab Emirates, April 2014

Middle East respiratory syndrome coronavirus (MERS-CoV) infection causes an acute respiratory illness and is associated with a high case fatality rate; however, the pathogenesis of severe and fatal MERS-CoV infection is unknown. We describe the histopathologic, immunohistochemical, and ultrastructural findings from the first autopsy performed on a fatal case of MERS-CoV in the world, which was related to a hospital outbreak in the United Arab Emirates in April 2014. The main histopathologic finding in the lungs was diffuse alveolar damage. Evidence of chronic disease, including severe peripheral vascular disease, patchy cardiac fibrosis, and hepatic steatosis, was noted in the other organs. Double staining immunoassays that used anti–MERS-CoV antibodies paired with immunohistochemistry for cytokeratin and surfactant identified pneumocytes and epithelial syncytial cells as important targets of MERS-CoV antigen; double immunostaining with dipeptidyl peptidase 4 showed colocalization in scattered pneumocytes and syncytial cells. No evidence of extrapulmonary MERS-CoV antigens were detected, including the kidney. These results provide critical insights into the pathogenesis of MERS-CoV in humans.

Two new Trechispora species from La Réunion Island

Two new species of corticioid Basidiomycota are described: Trechispora cyatheae Ordynets, Langer & K.H. Larss. sp. nov., and T. echinocristallina Ordynets, Langer & K.H. Larss. sp. nov. Both species were repeatedly collected in 2013 and 2015 in tropical rain forests of La Reunion Island (West Indian Ocean). In addition to La Reunion, we report T. echinocristallina from Papua New Guinea. Trechispora cyatheae has very small basidiospores not exceeding 3.5 μm in length, and sporulates on dead frond stipes of Cyathea glauca Bory, the tree fern endemic to La Reunion. Trechispora echinocristallina has thick dimitic fruit bodies found on tiny woody and non-woody forest debris and soil. It has unique crystal aggregations composed of thick-walled crystalline spheres, covered by numerous radially protruding needle-like crystals. This crystal configuration was hitherto not known in Trechisporales K.H. Larss. Both new species are illustrated as line drawings and microphotographs. Their relationships to other Trechispora species is evaluated by maximum likelihood and Bayesian inference of phylogeny based on concatenated ITS and LSU rDNA datasets.

Towards a phylogenetic generic classification of Thelypteridaceae: Additional sampling suggests alterations of neotropical taxa and further study of paleotropical genera.

Thelypteridaceae is one of the largest fern families, having about 950 species and a cosmopolitan distribution but with most species occurring in tropical and subtropical regions. Its generic classification remains controversial, with different authors recognizing from one up to 32 genera. Phylogenetic relationships within the family have not been exhaustively studied, but previous studies have confirmed the monophyly of the lineage. Thus far, sampling has been inadequate for establishing a robust hypothesis of infrafamilial relationships within the family. In order to understand phylogenetic relationships within Thelypteridaceae and thus to improve generic reclassification, we expand the molecular sampling, including new samples of Old World taxa and, especially, many additional neotropical representatives. We also explore the monophyly of exclusively or mostly neotropical genera Amauropelta, Goniopteris, Meniscium, and Steiropteris. Our sampling includes 68 taxa and 134 newly generated sequences from two plastid genomic regions (rps4-trnS and trnL-trnF), plus 73 rps4 and 72 trnL-trnF sequences from GenBank. These data resulted in a concatenated matrix of 1980 molecular characters for 149 taxa. The combined data set was analyzed using maximum parsimony and bayesian inference of phylogeny. Our results are consistent with the general topological structure found in previous studies, including two main lineages within the family: phegopteroid and thelypteroid. The thelypteroid lineage comprises two clades; one of these included the segregates Metathelypteris, Coryphopteris, and Amauropelta (including part of Parathelypteris), whereas the other comprises all segregates of Cyclosorus s.l., such as Goniopteris, Meniscium, and Steiropteris (including Thelypteris polypodioides, previously incertae sedis). The three mainly neotropical segregates were found to be monophyletic but nested in a broadly defined Cyclosorus. The fourth mainly neotropical segregate, Amauropelta, was found to include species considered to be part of Parathelypteris. In Old World thelypteroids, which correspond to nearly half the diversity in the family, an increase in sampling is still needed to resolve relationships and circumscription of genera, particularly in the christelloid clade (i.e., Amphineuron, Chingia, Christella, Pneumatopteris, Pronephrium, and Sphaerostephanos). Based on currently available knowledge, we propose the recognition of 16 genera in the family.